MADPH-08-1525 



Icosahedral (^5) Family Symmetry and the Golden Ratio 
Prediction for Solar Neutrino Mixing 

Lisa L. Everett and Alexander J. Stuart 
Department of Physics, University of Wisconsin, Madison, WI 53706, USA 

Abstract 

We investigate the possibility of using icosahedral symmetry as a family symmetry group in 
the lepton sector. The rotational icosahedral group, which is isomorphic to A5, the alternating 
group of five elements, provides a natural context in which to explore (among other possibilities) 
the intriguing hypothesis that the solar neutrino mixing angle is governed by the golden ratio, 
<ft = (1 + \/5)/2. We present a basic toolbox for model-building using icosahedral symmetry, 
including explicit representation matrices and tensor product rules. As a simple application, we 
construct a minimal model at tree level in which the solar angle is related to the golden ratio, 
the atmospheric angle is maximal, and the reactor angle vanishes to leading order. The approach 
provides a rich setting in which to investigate the flavor puzzle of the Standard Model. 



1 



INTRODUCTION 



The discovery [l| and subsequent measurements , |g, lZJ, ISj of neutrino oscillations 

have verified that neutrinos are massive and lepton mixing is observable, the first particle 
physics evidence of physics beyond the Standard Model (SM). The data have revealed an 
intriguing pattern of lepton flavor mixing parameters, providing a new facet to the SM flavor 
puzzle. Current measurements do not allow for an unambiguous extraction of the neutrino 
mass pattern. However, global analyses of the data [9j have provided the following 3a limits 
on the mass-squared differences Am-; = m 2 — mf {rrii are the neutrino mass eigenstates): 



Am 2 12 = 7.671 



:°f x 10- 5 eV 2 , 



(1) 



and 



Am 2 13 



(2) 



' -2.37^ x l(r 3 eV 2 (inverted hierarchy) 
. 2.46lg;42 x lCT 3 eV 2 (normal hierarchy) . 
Therefore, not only is there a suppression of the neutrino mass scale with respect to the 
electroweak scale, which is known both from direct bounds [lo| and cosmological constraints 



111] , but the neutrino oscillation data also suggest a relatively mild neutrino mass hierarchy 
as compared to the mass hierarchies of the charged fermions. 

The mixing angles of the Maki-Nakagawa-Sakata-Pontecorvo (MNSP) 



12 



131 ] lepton 



mixing matrix are also constrained to the following 3a ranges (in degrees): 







12 



42.3+^ 7 3 , 







13 



0.0 



+12.9 

-o.o • 



(3) 



In contrast to the small angles of the Cabibbo-Kobayashi-Maskawa (CKM) 14j quark mixing 
matrix, the MNSP matrix has two large angles, the solar angle (#12) arid the atmospheric 
angle (#23); and only one small angle, the reactor mixing angle #13, which is bounded to be 
less than the quark Cabibbo angle, the largest angle of the quark sector. 

To understand the origin of fermion masses and mixings in the SM, the paradigm is to use 
family symmetries to restrict the form of the mass matrices, and to generate mass hierarchies 
and mixings through small symmetry breaking perturbations. For the lepton sector, it is 



standard to work within the the framework of the seesaw mechanism [15|, which provides 
arguably the most natural explanation for the suppressed neutrino mass scale. In terms of 
mixing angles, it is natural to perturb about vanishing mixings for the quarks. However, the 
presence of large angles in the lepton sector hints at a range of possible theoretical starting 



points. The data suggest that at leading order, 9±3 is likely to be zero and #23 is likely to be 
maximal (45°), but there are a range of possibilities for 812. The approach known as "quark- 



lepton complementarity" [16|, which has been widely studied in the literature, assumes that 
Q\2 is also maximal at leading order, but is shifted to values within the experimental range 
by Cabibbo-sized effects. A perhaps even more popular approach is to seek the Harrison- 



Perkins-Scott (HPS) "tri-bimaximal" mixing pattern [17j at leading order, in which 

12 = ^2 PS) = tan" 1 (^=j = 35.26°, (4) 

well within the allowed range of the data. An attractive feature of the tri-bimaximal mixing 
scenario is that it emerges quite naturally from flavor theories based on non-Abelian discrete 
family symmetries which are subgroups of 5*0(3) or SU(3). Models have been constructed 
based on A4 (which is isomorphic to the tetrahedral symmetry group T) jl8j ]. the bin ary 



tetrahedral group (2) T [19], A(3n 2 ) [20j, the semidirect product of Z 3 and Z 7 [2l|, <S 4 [22] 
(which is isomorphic to the rotational symmetry of the cube), and S3 2s| (here additional 
input is needed to achieve the HPS value of #12), among others. The plethora of viable flavor 



models within this approach (see [24| for reviews) suggests that discrete non-Abelian family 
symmetries may be the key to understanding the flavor patterns of the lepton sector. 

In this paper, we explore the possibility that the rotational icosahedral symmetry group 
X, which is isomorphic to A5, the alternating group of five elements, is the family symmetry 
group in the lepton sector. Icosahedral symmetry is the only symmetry based on Plato's five 
"perfect" solids (the tetrahedron, cube, octahedron, dodecahedron, and icosahedron) which 
has not been used as a family symmetry group (the cube and the octahedron are described 
by the same symmetry group, as are the dodecahedron and icosahedron). This group has 
remained comparatively unexplored for physics applications, at least in part because it is 
not a crystallographic point group. Icosahedral symmetry is a logical candidate for a family 
symmetry, in that it has irreducible triplet representations (as do the groups given above, 
except for S3), and it contains the tetrahedral (A4) group as a subgroup. 



An intriguing feature of the icosahedral 



setting in which to explore the idea 25], l26|] that the solar mixing angle is related at leading 



'^.5) symmetry group is that it provides a natural 



order to the golden ratio, 4> = (1 + y/E)/2 (a solution of x 2 — x — 1 = 0): 

012 = tan" 1 (]■) = 31.72°, (5) 



which is 2a below the current best fit value (but can easily be shifted to higher values by 
small perturbations). This hypothesis was first mentioned in [2^], and recently explored in 



the context of Z 2 x Z 2 models in 



261 ] (it is worth noting that an alternative pattern in which 



cos 6*12 = <p/2 has also been suggested (27|). It was suggested in [26J that A$ might naturally 
generate this pattern, as the golden ratio is related to the geometry of the icosahedron. 

In this paper, we investigate the rotational symmetry group of the icosahedron as a family 
symmetry group, with the specific goal of building flavor models which result in the golden 
ratio prediction for the solar mixing angle as in Eq. (jSJ). To this end, we begin by reviewing 
basic properties of the rotational icosahedral symmetry group. We will identify a particularly 
useful group representation for our purposes, first given in the literature by Shirai 28(, and 
use it to construct tensor products of the irreducible representations. Finally, we apply this 
formalism by constructing toy flavor models which result in Eq. (jSJ) at leading order. Our 
studies will demonstrate that the icosahedral (^5) symmetry group provides a rich setting 
in which to explore the flavor puzzle of the SM. 



THEORETICAL BACKGROUND 



In this section, we will provide a brief overview of the basic properties of the icosahedral 
symmetry group. The methodology and many of the basic group p roperties are well known 



in the literature and can be found in 



28 



29, S, 



31 



32 



33 



34 



35, 



361 ] . However, the 



explicit tensor product decomposition has not to our knowledge been presented in the 



28| ) that we will 



literature, at least not for the specific group representation (due to Shirai 
argue is particularly well suited for the purposes of flavor model building. We compute 
these quantities at the end of this section. 



Preliminaries. The rotational icosahedral group, X, is the group of all rotations of the 
icosahedron which take vertices to vertices. Recalling that the icosahedron is the Platonic 
solid which consists of 20 equilateral triangles, such that it has 20 faces, 30 edges, and 12 
vertices, it is straightforward to see that there are five types of such rotations (we consider 
only proper rotations, and therefore ignore inversions). These rotations include rotations 
by or 2tt (the identity), rotations by 2ir/5 and Air/ 5 about an axis through each vertex, 
rotations by 27r/3 about axes through the center of each face, and rotations by it about the 



4 



midpoint of each edge. These five types of rotations form five disjoint conjugacy classes, 
which are denoted as follows: 



1, 12C 5 , 12C 5 2 , 20C 3 , and!5C 2 . 



(6) 



In the above, we follow the standard practice of using Schoenflies notation, in which C\ 
is a rotation by 2/c7r/ra, and the number in front identifies the number of group elements 
in the conjugacy class. It is worth noting here that this implies that elements of C\ are 
order n elements (elements which result in the identity after n operations). As can be seen 
from this classification, the rotational icosahedral group has 60 elements (it is order 60). 
Given that the order of the group is equal to the sum of the squares of the dimensions of the 
irreducible representations, and furthermore that the number of conjugacy classes is equal to 
the number of irreducible representations, the dimensions of the irreducible representations 
are given by 

1 + 12 + 12 + 15 + 20 = 60 = l 2 + 3 2 + 3 2 + 4 2 + 5 2 . (7) 

Hence, X has five irreducible representations: 1, 3, 3', 4, and 5 (note the presence of two 
distinct triplet representations, 3 and 3'). The single one-dimensional representation, as 
opposed to the case of Aj, which has three one-dimensional irreducible re pres entations, 

32l. 



results because X is a perfect group (it is equal to its commutator subgroup) [31 



The character table of the rotational icosahedral group is given in Table I (this table can 



also be found in 



28 



29 



31. 



33 



34 



35 



36j). The character table labels the character (trace) 



of the matrices which represent the group elements of a particular conjugacy class within a 
given irreducible representation. Table I indicates the appearance of the golden ratio for the 
triplet representations of the 12C 5 and 12C 2 conjugacy classes. Using the character table, 
the decomposition of tensor products of various irreducible representations can be deduced. 
More explicitly, the product of two three-dimensional representations 3©3=1©3©5 (and 
similarly for 3' ® 3', with 3 — > 3'). The full list of tensor products is presented in Table II (it 



also appeared earlier in the literature 



35]). Given the properties of tensor products of 5*0(3), 



the tensor product decomposition of two elements of the same irreducible representation can 
be classified by symmetry properties. For example, 3©3 = 1©3©5 = (1©5)s©3 j 4, in which 
the labels correspond to symmetric or antisymmetric, respectively (see {35] for the complete 
listing). We also show in Table III the relations between the irreducible representations of 



X and the irreducible representations of A$, which are also found in 



34]. 



ill which (p = (1 + y/E)/2. 
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1 


3 


3' 


4 


5 
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1 


3 


3 


4 


5 


12C 5 


1 




1-0 


-1 





YlCl 


1 


1-0 


4> 


-1 





20C 3 


1 
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-1 


15C 2 


1 


-1 


-1 
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: the icosahedral symmetry group 
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3' 
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3' =1 
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4© 5 
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<8 


)4=3' 


©4© 


5 
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'0 


S>4=3 


©4© 


5 








30 5 


=3© 


3' ©4 


9 5 






3' (5 


5 5=3 


3' ©4 


© 5 






4(8)4= 


1©3( 


e3'© 


49 


35 






4©5 




3 3' 


©4© 


5 9 


35 




5$ 


5 5=1 


9 


3 3' 


©4© 


49 


35 


© 5 



TABLE II: The tensor product decomposition for the rotational icosahedral symmetry group X 
and the alternating group A5; see also [351 ]. 



Group Presentation and the Shirai Basis. As is the case for any finite group, the 
elements of the rotational icosahedral symmetry group can be generated by a set of basic 
elements which satisfy certain relations; the elements and rules together are known as the 
"presentation" of the group. Unlike the case for smaller groups (such as A4), there are 
several equivalent ways to present the group, which either involve two generators or three 
generators; a classification can be found in [29Q . We will proceed here by choosing a specific 
presentation of X due to Hamilton (see e.g. 

<a,b\a 2 = b 3 = {ab) 5 = e> . (8) 





A± 


1 


1 


3 


3 


3' 


3 


4 


301 


5 


3 © 1' © 1" 



TABLE III: The decomposition of the irreducible representations of the rotational icosahedral 

n 

symmetry group X in terms of the irreducible representations of A4; see also [34]. 



This notation indicates that with the two generators, a and b, the entire rotational icosa- 
hedral group can be generated (or something isomorphic to X) when a and b are combined 
subject to the "rules" of Eq. (jSJ). Here e is the identity, as is standard in the mathematical 
literature. An explicit example of this presentation is given in 35 1. 

The two generators of Eq. (jHJ) belong to the C2 and C3 conjugacy classes, as they square 
and cube to the identity, respectively. We have found that for the purposes of flavor model 
building with the golden ratio prediction for the solar angle in mind, it is advantageous to 
use a presentation which explicitly uses an element either from C5 or Cf, as the golden ratio 
appears in the character table for these conjugacy classes explicitly (see Table I). As noted 



e.g. m 



361 ]. Hamilton's presentation can be rewritten in terms of an order two generator 



S = a and an order five generator T = bab, such that 



< S, T\S 2 = T 5 = (T 2 ST 3 ST' 1 ST ST' 1 ) 3 = e> . 



(9) 



xplicit representations were presented both by Hoyle 36) and Shirai 28j. Hoyle's basis 



361 ] is not unitary for a subset of the irreducible representations, and hence is not usefu 
flavor model building. In Shirai's basis, which incorporates generators from C5 and C2 
the generators for the triplet representation (3) follows: 



for 



28|, 



/ 



-1 



1 \ 



i 1 



1 



1 \ 



-d> 4 1 



(10) 



7 



while for the other independent triplet representation (3'), the generators take the form 



5, 



V 



T,, 



( 



\ 



\ 



(11) 



The generators for the four- dimensional representation in the Shirai basis are given by 



\ 



( 



Ta 



J 



\ 



-1 


1 


-3 


75 


-1 


-3 


1 


75 


3 


1 


1 


75 


75 


-75 


-V5 


-1 



( -1 -1 -3 -75 \ 
-1 3 1 -75 
-3 1 -1 75 
-75 -75 75 -1 

while for the five- dimensional irreducible representation, the order two generator is 

s/3 \ 

40 

yj 

20 

VM 

2 
2 



(12) 



4 

ii 

2 

1 

202" 

yj 

2 

\ 40 



2 

2 


1 

20 2 


V5 
2 


1 


1 





1 





-1 





-1 


1 


V3 
20 


V3<f> 
2 


V3 
2 



(13) 



and the order five generator is 



/ 1^30 
4 

iL 

2 
1 

202" 
2 

vj 

40 



2 

-1 
1 



20 



l V5 yj \ 

" 202 2 40 

yj 

20 

yj0 

2 



1 

1 







1 



1 yJ 

X 2 



yj30 30-1 
2 2 4 



(14) 



(our Eqs. (fT3l) - (fT4l) correct minor typos found in [28|). Before proceeding further, let us 
first comment on other presentations that are available in the literature. The icosahedral 
(.4.5) group can be defined using three generators, in which two of the generators define the 
tetrahedral (^.4) subgroup. Such a presentation is given in the classic text of Lomont [31[ 
and the work of Cummins and Patera [34], as follows: 



< A 1 ,A 2 ,A 3 \A 3 1 



A\ 



A 2 3 = (A 1 A 2 f 



(A 2 A 3 ) 3 



e > 



(15) 



In the above, the generators A\ and A 2 generate the tetrahedral subgroup of X. The explicit 
representations, which we will not state explicitly here, differ in 31] and (3^] by a similarity 



8 



transformation. In terms of the generators of Hamilton's presentation a and b, we have the 
following equivalence relations (up to a similarity transformation): 

A 1 ~ [b(ab 2 ) 2 ab(ab 2 ) 2 } 2 , A 2 ~ a, A 3 ~ 6a6 2 a[6(a6 2 ) 2 a] 2 [6(a6 2 ) 2 ] 2 . (16) 

While any of these presentations can be used as a basis for flavor model building, the Shirai 
basis has the advantage that the golden ratio appears explicitly, which will turn out to be a 
useful feature for seeking the golden ratio prediction for the solar mixing angle as given in 
Eq. (jSJ). Therefore, in this paper we will use the Shirai basis (Eqs. ffTOj) — ffT^jl ) . The explicit 
tensor product decomposition has not, to our knowledge, been presented in the literature 
for this basis. We construct it explicitly in the next section. 



Tensor products. The next task is to find group invariants of X. General procedures for 



constructing the invariants for finite groups can be found in |30l.l3l|. with explicit applications 
to icosahedral (.4.5 ) symmetry in 34],|35j]. In this paper, we will focus for the sake of simplicity 
on constructing the relevant Kronecker products within the Shirai basis. Our results are 
sufficient for flavor model building at tree level using this basis. Let us begin by defining 
the two distinct triplet states as 3 = (aj, a 2 , a 3 ) T and 3 = (b\, b 2 , b 3 ) T . Using Eqs. (TTOT l (TT^j) . 
it is straightforward to show that the decomposition of 3 ® 3 = 1 © 3 © 5 yields the following 
results for the singlet, 

1 = ai&i + a 2 b 2 + a 3 b 3 , (17) 

the triplet, 



3 = (a 3 b 2 — b 2 a 3 , a x b 3 — a 3 bi, a 2 b\ — ai& 2 ) 



(18) 



and the fiveplet: 



a 2 b 2 — ai&i 
a 2 bi + a x b 2 
a 3 b 2 + a 2 b 3 
aib 3 + a 3 b x 

\ ~73^ ai&1 + fl2&2 ~ 2fl 3 & 3) ] 

which show the symmetry of the singlet and the fiveplet, and the antisymmetry of the triplet. 
The Kronecker products for 3' <S> 3' for the singlet and the triplet have an identical structure 



9 



to those of the 3 <S> 3, but for the fiveplet we have the different symmetric combination: 



l (-4ai&i - 4>a 2 b 2 + V5a 3 b 3 



\ 



a 2 bi + ai& 2 
-(a 3 6i + a x b 3 ) 

\ ^-((1 - 30)aA + (30 - 2)a 2 b 2 + a 3 b 3 ) ) 
For 3 ® 3' = 4 © 5, we have the following result for the fourplet: 



\a 3 b 2 - 4>atb 3 



\ 



'i + ^a 2 b 3 
-\o-ibi + <M 2 6 2 
a 2 bi - a x b 2 + a 3 b 3 



(19) 



and the five- dimensional representation takes the form: 

' l((f) 2 a 2 bi + j?aib 2 - V5a 3 b 3 ) ^ 
-(0ai&a + ^a 2 &2) 
ia 3 6i - (j)a 2 b 3 ■ (20) 

V 2 (^ a 2&i + <\>&\b 2 + a 3 b 3 ) J 

For our purposes in this paper, the only remaining Kronecker products we need are (i) 
the singlet which results from tensoring together two fourplets, 4 = (a x , a 2 , a 3 , a 4 ) T and 
4 = (pi, b 2 , b 3 , b i ) T , which is given by 



1 = a x b x + a 2 &2 + &zb 3 + a 4 6 4 , 



(21) 



and (ii) the singlet which results from the tensor decomposition of the direct product of two 
fiveplets, 5 = (ai, a 2 , a 3 , a 4 , a 5 ) T and 5 = (pi, b 2 , b 3 , 64, & 5 ) T ', which takes the familiar form 



1 = a x b x + a 2 b 2 + a 3 6 3 + a 4 6 4 + a 5 6 5 . 



(22) 



The form of these sing lets results directly from unitarity and the reality of the Shirai basis. In 
forthcoming work 37(], we will present a more comprehensive analysis of the group invariants 
which connects with the general formalism presented in 34j, [35 ] . 
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3 3 


3 


3' 


3 


3' 


3' 


3' 
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3 3 


3' 


3 


3' 


3 


3' 


3' 


N 


3 3' 


3 


3 


3' 


3' 


3 


3' 



TABLE IV: Possible charge assignments under I for the SM lepton doublets L, lepton singlets e, 
and the right-handed neutrinos N. 

LEPTON FLAVOR MODEL BUILDING 



In this section, we apply this formalism to the issue of flavor model building based on the 
rotational icosahedral symmetry group as a family symmetry of the lepton sector. We focus 
in this paper exclusively on the lepton sector and defer the quark sector to future work. 
With this restriction in mind, the starting point is to examine possible family symmetry 
charge assignments for the SM fields, which include the lepton doublets L iy the charged 
lepton singlets e^, and, in the case of seesaw models, the right-handed neutrinos Ni (i is 
a family index). The electroweak Higgs field(s) are assumed to be blind to the family 
symmetry. Clearly, it is natural for the lepton sector fields to transform as a 3 or a 3' under 
the icosahedral symmetry. This leads to eight possible scenarios, as shown in Table IV. 

For concreteness, we will focus here on the case in which L is a 3 and e is a 3', leaving 
the A^ charge unspecifed for the moment. With these charge assignments, the standard 
procedure is to use the Froggatt-Nielsen mechanism 38( and introduce flavon fields, scalar 
fields charged under the family symmetry which acquire vacuum expectation values (vev's) 
that break the family symmetry. To determine the options for the flavon field content, recall 
that the structure of the mass terms we aim to generate takes the standard schematic form: 



(23) 



-C rn = Y^LiNjH + M^Nify + Y^LiZjH + h.c, 

in which Y^ denote effective Yukawa couplings, Ml- represents the Majorana mass matrix 
of the right-handed neutrinos, and H denotes the electroweak Higgs field. (For simplicity, 
we assume a single Higgs field, as in the SM; it is straightforward to generalize the analysis 
for the case of two electroweak Higgs doublets, as in the minimal supersymmetric extension 
of the SM.) The Dirac mass matrices Mjp are then given by Y^(H). Within the seesaw 
framework, the eigenvalues of are much larger than the Higgs vacuum expectation 

value, such that the heavy right-handed neutrinos can be integrated out to yield the standard 



11 



seesaw mass matrix for the light (primarily left-handed) neutrinos 15]: 



M v = M D (M (M) )" 1 M^. (24) 



This mass term can also be encoded by Weinberg's dimension five operator 39| : 



-C w = ^LiHLjH, (25) 

where M is a (presumably high) cutoff scale (i.e., the scale of the eigenvalues of the Major ana 
mass matrix of the right-handed neutrinos in a seesaw framework). For simplicity, we 

will use Eq. (1251) to encode the effective seesaw mass term, which bypasses the details of the 
right-handed neutrinos, for the time being. We will return to this issue later in the paper. 
In our scenario, the effective mass terms are thus of the LL and Le type, with 

LL : 3 ® 3 = 1 © 3 © 5, Le : 3 © 3' = 4 © 5. (26) 

We see that in the absence of flavor symmetry breaking, the bare LL term is allowed, but 
the bare Le term is absent. Given the structure of the singlet coupling of two triplets (see 
Eq. ( JT7j) ). this indicates that the neutrinos can have degenerate masses of order (H) 2 /M and 
the charged leptons are all massless at leading order. This result is, of course, quite far from 
the observed pattern of lepton flavor mixing. However, the flavor structure emerges through 
the spontaneous breaking of the icosahedral symmetry through the vacuum expectation 
values of flavon sector fields. This is a generic feature for flavor models based on non-Abelian 
discrete family symmetries. The details of the flavon sector matter content and couplings are 
thus the key elements to consider for flavor model building. Within icosahedral symmetry, 
we see that the flavon sector fields of interest are singlets and fiveplets for the neutrino 
sector, and fourplets and fiveplets in the charged lepton sector (the symmetry of the LL 
term disallows the coupling of LL to a flavon field which transforms as a triplet). 

Even within this simpler class of scenarios in which the L and e fields have been assigned 
to specific representations, to address whether the golden ratio prediction for the solar 
mixing angle can naturally emerge within such icosahedral (A5) models requires a systematic 
study, which is currently in progress 37|]. Here, for illustrative purposes, we will focus on 
constructing a minimal (toy) scenario which reproduces a similar structure of mixings as 



found in 



261 ]. in which the solar mixing angle originates from the neutrino sector, and the 



maximal atmospheric mixing angle arises from the charged lepton sector. In this paper, we 

12 



do not attempt to construct a fully realistic model; instead, the goal is to see, through this 
very simple example, what flavon sector fields and associated vacuum expectation values 
are needed to reproduce this pattern within the framework of icosahedral symmetry. To 
this end, we include a "minimal" number of flavon fields and only consider leading order 
couplings. Furthermore, we will not include effects of renormalization group running in this 
scenario, which would be necessary in a comprehensive study (see e.g. 401]). 

The simple model we consider has the following flavon sector fields: the fields £ and i[>, 
which each transform as a 5, and the flavon Xi which transforms as a 4. The presence of 
two distinct fiveplets turns out to be unavoidable in this context; we will comment on this 
shortly. We will assume that £ couples to the LL operator and that x an d ^ couple to the 
Le operator; this assumption necessitates additional symmetries which forbid £ and ip from 
coupling to Le and LL, respectively. For simplicity, we will also assume that the tree level 
LL term is forbidden by such additional symmetries. This tree level term is proportional to 
the identity matrix, and hence it does not affect the mixing angles or the predictions for the 
mass-squared differences, but does affect the overall neutrino mass scale. Similarly, singlet 
flavon fields which couple to LL would also just affect the overall neutrino mass scale. We 
assume such singlet fields are absent. 

With these assumptions, the mass terms take the following schematic form: 



r 



fiijk 



lijl 



-LiHLjH^k + -^-LiejH^k + -j-L^Hxi + h.c 



(27) 



MM' J ™ M' M' 
in which M' represents the (presumably high) scale of flavor symmetry breaking, and a^, 

Pijk, and jiji are dimensionless couplings that encode the tensor product decomposition of 

the icosahedral symmetry as well as any (presumably 0(1)) coupling factors. In principle, 

M' has no relation to M, particularly within the context of the neutrino seesaw. After 

electroweak symmetry breaking, the resulting mass matrices are given by 

£2 £1 — £3 



M„ 



a 



V 



(28) 



6 75^5 / 



Me 



-H2 ^1 + ^f>5 



^1 + f ^5 



^2 



^3 



+ 7 



4xs -X4 -</>xi \ 



X4 



5X1 



0X2 

X4 J 



(29) 
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in which a, (3, and 7 are dimensionless couplings (functions of the parameters of Eq. ( 1271) ) 
which encode the suppression of the neutrino and the charged lepton mass scales. 

To determine the nature of the flavor symmetry breaking and the resulting masses and 
mixings, the dynamics of the flavon sector must be considered in detail. The flavon structure 
is clearly quite rich due to the presence of the four- and five-dimensional representations. 
Given the complexity of this sector, we will leave the flavon dynamics unspecified in this 
paper as a first step toward flavor model building, and defer this important issue for future 
work. We will instead examine the possible vev patterns that can result in scenarios anal- 
ogous to those presented in 26], and defer a more systematic study of other viable options 
to future work. Within this approach, the golden prediction for the solar neutrino mixing 
angle can be obtained in the neutrino sector by assigning the following vev to £ : 







' vm( m2 ~ mi ) N 




— ► ^ 


^(roa-mi) 
















Us/ 




^ -(mi + m 2 ) / 



(30) 



In the above, m\ and m 2 are mass parameters, which turn out to be the neutrino mass 
eigenvalues associated with the solar neutrino pair. The dominant component of the vacuum 
expectation value of £ is that of £5; the other terms represent small corrections, of order 
~ Amf 2 /('7ii + From Eq. ( l28l) . one sees that the vanishing of the £3,4 component vevs 
is needed to ensure that only the solar mixing angle results from the neutrino sector. The 
solar angle in this scenario is determined by the interplay between the £1 and £ 2 vevs, as the 
£ 5 vev leads to a diagonal mass matrix. Explicitly, the neutrino mass matrix is 

' 0mi + ^m 2 m 2 — mi ^ 

m 2 — mi imi + 0m 2 • (31) 

\ — \/%(mi + m 2 ) j 

The neutrino mass eigenvalues are mi, m 2 , and m 3 = —(mi + m 2 ). The constraint on 
results from the interplay of the small £1 and £ 2 vevs with the £5 vev, which is needed to 
obtain the golden ratio prediction for the solar mixing angle. The negative sign of m 3 can be 
absorbed by a phase rotation because the neutrinos are Majorana fermions. The neutrino 



My 
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mixing matrix U v , defined by U v M v Ul ', satisfies Eq. 



by design: 



1 

V5<t> 



^ o \ 

I A. n 

V5 U 



V 



(32) 



/ 



-z 

The masses mi and m2 are determined by accommodating the splittings of the solar and 
atmospheric neutrino mass-squares (see Eqs. ([I])-©); we obtain mi = 2.77 x 10~ 2 eV, 
m 2 = 2.91 x 10 -2 eV, and |m 3 | = 5.68 x 10~ 2 eV. Hence, the resulting pattern has a 
normal mass ordering, but with a mass spectrum which is quasi-degenerate. The mass pattern 
obtained here differs from that of |26J], in which was an independent parameter, while 
mi and m 2 were given by m 1 = m/<fi, m 2 = 4>m (Am 2 2 = y/Em 2 in their scenario). 

In the charged lepton sector, we seek maximal mixing between the second and third 
families. At leading order in flavon fields, we assume only m T is nonvanishing, and that the 
charged lepton mass matrix takes the form 

/ N 



1 

V2 



as m 



m T 
V m T / 



(33) 



26j. Such a pattern can be achieved with the following vev assignments for \ an d ip- 



^3 

/ 

Xi 

X2 
X3 



2^/3 



3 




1 



(34) 



3^7 



(35) 



V 1 / 



Given Eq. (I3B1 . the vevs of ?/> 2 ,4 and Xi,3 are required to vanish due to the structure of 
Eq. (|29l) . The left-handed states are diagonalized by the mixing matrix 



\ 



1 
-i L 



V° ts is J 



(36) 
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It can be shown that Eq. (|33|) cannot be reproduced with nonvanishing vevs for only \ or 
only ip (see Eq. ( 129]) ). Furthermore, two distinct fiveplets are required in this model (one 
which couples to the neutrino sector and one which couples to the charged lepton sector), 
since the required vev of ip is qualitatively different than that of £ (see Eq. (|30p ). 

Nonvanishing electron and muon masses can also be accommodated by fitting the full 
result for the charged lepton mass matrix given in 26j: 

( V2m e ^ 





1 

V2 



V 



/ 



by allowing for nonvanishing vevs for -02,4 and Xi,3 as follows: 



^4 



^(-(j)y/2m e 



2 4>m r 



(37) 



(38) 



X2 
X3 

\ XA ) 



7 



V 



-\m, 



m e + 0m p 



/ 



(39) 



Such small entries in the flavon vevs clearly represent a large fine-tuning, which is needed to 
achieve the hierarchically suppressed electron and muon masses at leading order in flavon 
fields. In general, whether or not flavon vevs such the ones presented above can emerge nat- 
urally within icosahedral symmetry is a question which can only be addressed by explicitly 
considering the dynamics of the flavon sector. Though the vacuum alignment issue must be 
addressed, it is not implausible that the appearance of the golden ratio in the flavon vevs 
can emerge naturally from the underlying icosahedral symmetry. 

We pause here to note that generating charged lepton masses and mixing is quite different 
here than in the case of Aa models, where typically the charged lepton singlets are assigned 
to the three different one-dimensional irreducible representations. In such A\ models, the 
hierarchical charged lepton masses can be addressed via separate small Yukawa couplings, 
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whereas in the X (A5) scenarios studied here, one necessarily has to address all three families 
together. In forthcoming work, we will analyze details of higher-order couplings which can 
result in the charged lepton mass matrix of Eq. (1371) . among other possibilities. 

The last feature remaining is to present the explicit form of the MNSP matrix. Given 
Eq. (|H2D and Eq. (J3SD, the MNSP matrix takes the form 



U = Umnsp = U e Ut 



( ,/x or o \ 

V V5 V V54> U 

1 r~r~ 1 rv i_ 

V2\ VE<t> V2V V5 V2 



V, (40) 



where V = Diag(l, 1, i) is the Majorana phase matrix. By design, the MNSP matrix has 
a vanishing reactor mixing angle, a maximal atmospheric mixing angle, and a solar angle 
given by Qyi = tan -1 (1/0), as in Eq. (jSJ). 

In this scenario, once the charges of L and e are specified and the flavon field content is 
chosen, there are in principle fourteen complex parameters (the flavon vevs) which are inputs 
for the model. As there are eight physical observables (3 charged lepton masses, two neutrino 
mass-squared differences, and three lepton mixing angles), without further input the model 
is under const rained. With the choice of flavon field vevs of Eqs. ( 1301) . (134"1) . and (1351) . this 
number of parameters has been reduced by hand to 8 real parameters. Therefore, we can 
predict observables such as the neutrinoless double beta decay parameter, mpp = J2i m iUei, 
which is given by 

mi(f) m 2 

mm = ^ + wf; (41) 

In our model, mpp = 2.8 x 10 -2 eV. In regards to CP phases, the Dirac MNSP phase cannot 
be predicted because #13 = to this order of approximation. As seen in Eq. (liOjl . our 
scenario does have a Majorana phase, but it does not enter into the expression for mpp. 

We close this section by commenting on seesaw realizations of scenarios such as the one 
presented above. Given the form of Eq. (1231 . we see that whether the right-handed neutrinos 
are assigned to transform as a 3 or a 3', the possible flavon fields in the right-handed neutrino 
sector are either a 1 or a 5 due to the symmetry of the NN operator. One straightforward 
example for the case in which is a 3 (more precisely, when L and N transform identically 
under X) is that the only flavon field which couples to LN and NN is the fiveplet field £. 
We will also assume that the tree-level coupling of NN is forbidden (e.g. by an additional 
symmetry). In this case, is proportional to the Dirac mass matrix, Mo, which results 
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from the LN coupling. The seesaw formula then results in 

M v = M d (M (m) ) _1 M£ ~ MdM^MI ~ M D , (42) 

since Md is a symmetric matrix in this case. If £ has the same vev pattern as in Eq. (|30|) . it is 
straightforward to see that we obtain exactly the same flavor mixing structure of the neutrino 
mass matrix of Eq. (l3Tj) . Therefore, this scenario represents a simple seesaw extension of the 
effective theories given above. However, it is clear that there are many other possible seesaw 
embeddings that are worthy of further exploration. We plan to analyze this and other crucial 
model-building issues, most notably the dynamics of the flavon sector, in future work. 



CONCLUSIONS 



In this paper, we have investigated the possibility that the rotational icosahedral sym- 
metry group X, which is isomorphic to As, is the family symmetry group of the lepton 
sector. Icosahedral symmetry has remained relatively unexplored for physics applications, 
most notably because it is not compatible with translational invariance and thus is not a 
crystallographic point group. Nevertheless, the rotational icosahedral symmetry group has 
many features which make it an appealing candidate for a family symmetry. In particular, 
the presence of two distinct triplet irreducible representations provides a rich framework in 
which to generate the large mixing angles of the lepton sector. 

We have focused in this paper on developing the tools to construct models in which 
the solar mixing angle is related in a nontrivial way to the golden ratio, = (1 + y/E)/2. 
Icosahedral symmetry provides a natural setting in which to explore this hypothesis, as 
the golden ratio emerges from the geometrical properties of the icosahedron. Although the 
golden ratio prediction for the solar mixing angle is an intriguing option, it is not the only 
plausible outcome in models with icosahedral symmetry. One possible research direction is 
to study how tri-bimaximal mixing might arise within this framework. This outcome is also 
in some sense natural, in that the tetrahedral (Aa) group is a subgroup of I. 

Our work represents the first stage for flavor model building using icosahedral symmetry. 
We have constructed explicit tensor products using a specific group representation due to 
Shirai [28|, and used these results to construct toy models which reproduce the golden ratio 



prediction for the solar mixing angle of Eq. (jSJ). Much work still needs to be done within this 
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framework, including a more general analysis of higher-order invariants and the implications 
for flavor model building beyond leading order, particularly with respect to the important 
issue of the flavon sector dynamics. Though our investigation is still at a preliminary stage, 
icosahedral symmetry clearly provides a rich arena in which to explore the flavor puzzle of 
the Standard Model and is certainly worthy of further study in this exciting era of continued 
experimental probes of lepton flavor mixing. 
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